Crystallization kinetics of nanoconfined GeTe slabs in GeTe/TiTe2-like superlattices for phase change memories

Superlattices made of alternating blocks of the phase change compound Sb2Te3 and of TiTe2 confining layers have been recently proposed for applications in neuromorphic devices. The Sb2Te3/TiTe2 heterostructure allows for a better control of multiple intermediate resistance states and for a lower drift with time of the electrical resistance of the amorphous phase. However, Sb2Te3 suffers from a low data retention due to a low crystallization temperature Tx. Substituting Sb2Te3 with a phase change compound with a higher Tx, such as GeTe, seems an interesting option in this respect. Nanoconfinement might, however, alters the crystallization kinetics with respect to the bulk. In this work, we investigated the crystallization process of GeTe nanoconfined in geometries mimicking GeTe/TiTe2 superlattices by means of molecular dynamics simulations with a machine learning potential. The simulations reveal that nanoconfinement induces a mild reduction in the crystal growth velocities which would not hinder the application of GeTe/TiTe2 heterostructures in neuromorphic devices with superior data retention.

that the nanoconfinement of Sb 2 Te 3 reduces the structural relaxations responsible for the resistance drift in the amorphous phase 7 .
However, pure Sb 2 Te 3 suffers from an inadequate data retention due to insufficient stability of the amorphous phase.It would therefore be of interest to investigate the possibility to substitute Sb 2 Te 3 with a phase change material with a higher crystallization temperature such as GST or GeTe.In this respect, we must however consider that nanoconfinement could have several possible effects on the crystallization kinetics 8 .In the case of elemental Sb, for instance, it was shown that the amorphous phase can be dramatically stabilized against recrystallization in ultrathin films 3-10 nm thick capped by insulating layers 9,10 .It is therefore important to assess whether the materials could retain a high crystallization speed at high temperature along with good data retention at low temperature when it is confined in ultrathin films.
On these premises, in this article we report on molecular dynamics (MD) simulations of the crystallization of the phase change compound GeTe in a nanoconfined geometry.The simulations aimed at addressing the effect of confinement of GeTe slabs in a superlattice made of alternating layers of GeTe and TiTe 2 similarly to the Sb 2 Te 3 /TiTe 2 superlattices of Ref. 7 .We have chosen GeTe because it shares many properties with the most commonly studied GST alloy and because of the availability of an interatomic potential for large scale simulations that we have previously devised in Ref. 11,12 .This potential was generated by fitting with a Neural Network (NN) method 13 a large database of energies of small models computed within Density Functional Theory (DFT).We have not explicitly considered a real TiTe 2 slab which would require to device a proper NN potential, but we mimicked the confinement effects of TiTe 2 on GeTe by a suitable capping as discussed in the next section.The simulations reveal that the crystal growth velocity in the confined geometry is indeed lower with respect to the bulk, but to an extent (about a factor two at the temperature of maximum speed) which does not hinder the exploitation of GeTe/TiTe 2 superlattice for neumorphic applications.

Results and discussion
At normal conditions, GeTe crystallizes in a trigonal geometry (space group R3m) 14 .This phase named α-GeTe, with two atoms per unit cell, can be viewed as a distorted rocksalt geometry with an elongation of the cube diagonal along the [111] direction and an off-center displacement of the inner Te atom along the [111] direction giving rise to a 3+3 coordination of Ge with three short and stronger bonds (2.84 Å) and three long and weaker (3.17 Å) bonds.In the conventional hexagonal unit cell of the trigonal phase, the structure can be also seen as an arrangement of GeTe bilayers along the c direction with shorter intrabilayer bonds and longer interbilayers bonds.The interplanar distance within a bilayer is 1.506 Å while the interplanar distance across two bilayers is 2.062 Å.The trigonal ferroelectric phase transforms into the cubic paraelectric phase ( β-GeTe, space group Fm 3 m) above the Curie temperature of 705 K 15 .In the cubic phase, the alternation of long and short bonds survives in a disordered manner along all equivalent <111> directions as revealed by extended x-ray absorption fine structure (EXAFS), x-ray total diffraction measurements 16,17 and MD simulations 18 .However, more recent molecular dynamics simulations 19 suggest that the order-disorder character of the phase transition is weaker than as inferred from EXAFS data.
We have simulated the effect of confinement on the crystallization of GeTe by considering a slab made of nine bilayers of α-GeTe (18 atomic planes with a thickness of about 3 nm), encapsulated by capping layers aiming at mimicking the confining slabs of TiTe 2 in GeTe/TiTe 2 superlattices.The capping layer mimicking TiTe 2 on each side is made by a frozen bilayer of crystalline GeTe itself constrained at the lattice parameter of TiTe 2 as shown in Fig. 1a.This choice is made because for the MD simulations we used a NN potential suitable to describe only Ge-Ge, Ge-Te and Te-Te interactions.In fact, TiTe 2 is a layered hexagonal crystal (space group P 3m1) made of trilayer Te-Ti-Te blocks stacked along the c axis and separated by van der Waals gaps.The geometry of the hexagonal Te layers is the same in TiTe 2 and GeTe albeit with different lattice parameters, namely a = 3.7795 Å for TiTe 2 20 and a = 4.1677 Å for GeTe 14 .A good commensuration between a trilayer of TiTe 2 and α-GeTe in the hexagonal xy plane is obtained by considering multiples of the orthorhombic supercells with edge a and √ 3a , namely 11 x 10 orthorhombic cells of TiTe 2 and 10 x 9 orthorhombic cells of α-GeTe.The misfit is only 0.2 % along x and 0.7 % along y.We finally set the in-plane lattice parameters of the supercell to those of α-GeTe which means that the bilayers mimicking TiTe 2 are slightly strained by the amount given above.The model thus contains 10 x 9 x 2 = 180 atoms per atomic layer of GeTe and 11 x 10 x 2 = 220 atoms in each of the bilayers mimicking TiTe 2 , for a total amount of 4120 atoms.Periodic boundary conditions are applied along the three cartesian axis (see Fig. 1a).The TiTe 2 -like bilayers are oriented in such a way to expose the Te layer to the α-GeTe slabs on both sides as it would occur for a TiTe 2 slab.The distance between the two TiTe 2 -like bilayers is fixed to an arbitrary value of 4 Å as the two bilayers are just meant to mimic TiTe 2 slabs encapsulating GeTe from both sides and they must not interact with each other.The distance between the TiTe 2 -like bilayers and GeTe is instead optimized by minimizing the energy with respect to the c axis of the supercell.However, since this is just the initial model to be subjected to a thermal cycle of amorphization/crystallization, we decided to constraint the distance between GeTe and the Te layer of the capping layers to be the same on both sides, although on one side we have a Ge-Te contact and on the other side a Te-Te contact (see Fig. 1a).This choice is made because after a thermal cycle of amorphization/crystallization, we expect that the amorphous GeTe would face in the same manner the capping layer on both sides.We will discuss further this issue later on.This procedure finally yields an interplanar distance between Te atoms of the capping layers and the outermost atomic layer of GeTe of 3.463 Å which is close to interplanar distance of 3.48 Å between TiTe 2 and Sb 2 Te 3 in the TiTe 2 /Sb 2 Te 3 superlattice of Ref. 7 .We will refer to this model of GeTe encapsulated by TiTe 2 -like capping layer as the superlattice (SL) configuration.Because of the large lattice mismatch between GeTe and TiTe 2 , we assume that crystal nucleation would not be triggered by the interaction with the TiTe 2 slabs.Under this assumption our fake TiTe 2 would be sufficient to mimic the confinement effects introduced by TiTe 2 .As we will see in the following, crystal nucleation still occurs at the surface of the amorphous slab, not because of interaction with TiTe 2 but because of an atomic layering at the surface of the amorphous phase.
Since TiTe 2 has a much higher melting temperature than GeTe, we mimicked the confinement by TiTe 2 by freezing the atoms of the crystalline GeTe-like capping bilayers during the thermal cycle.Amorphization of the GeTe slab was achieved by equilibrating the system first at 1500 K for 200 ps and then at 1000 K for 100 ps.The liquid-like slab was then quenched to 300 K in 100 ps.
A snapshot of the amorphous GeTe (a-GeTe) slab encapsulated by the frozen capping bilayers mimicking TiTe 2 (superlattice configuration) is shown in Fig. 1b.We observed a small expansion of the amorphous slab which leads to a density of about 0.0363 atom/Å 3 to be compared with the initial density of 0.0372 atom/Å 3 of the slab in the crystalline phase.Because of the constraints imposed by the capping layers, the density of the amorphous slab is, however, sizeably larger than the theoretical equilibrium density of a-GeTe of 0.0315 atom/Å 3 given by the NN potential 21 .We mention that the theoretical density of a-GeTe computed with the NN potential (and consistently with the DFT framework) 21 is lower than the experimental value of 0.03327 atom/Å 322 .
We compared the structural properties of the amorphous phase to those of a bulk model quenched from the melt at the density fixed to that of the crystalline α-GeTe of 0.0372 atoms/Å 3 .Structural properties of the resulting amorphous slab were computed over 40 ps simulation at 300 K. Partial pair correlation functions, bond angle distribution functions, and distribution of the coordination numbers for the a-GeTe slab are compared with the results for the bulk in Fig. 2. The average partial coordination numbers of the slab and the bulk are compared in Table 1.The coordination numbers are obtained by integrating the partial pair correlation functions up to the bonding cutoff of 3.00 Å for Ge-Ge and Te-Te bonds and 3.22 Å for Ge-Te bonds.The coordination numbers are lower in the superlattice than in the bulk in part because of the slightly lower density (0.0363 atom/Å 3 vs.0.0372 atom/Å 3 ).However, there is a further reduction in the coordination number of Te because of an enrichment of Te at the two surfaces of the amorphous slab facing the Te planes of the capping layers.This feature is reproduced also by repeating the simulation for two different models with a shorter distance between the Te plane of the capping layer and the outermost plane of the GeTe block (see Fig. S1 in the Supplementary Information).In a first model, we changed the orientation of the TiTe 2 -like GeTe bilayers in such a way that the GeTe slab faces a Te capping layer on one side and a Ge capping layer on the other side both at a distance of 3.03 Å.In a second model, the interplanar distances between the Te layer of the TiTe 2 -like GeTe bilayer and the Ge and Te layers on the two side of the GeTe slab were obtained from DFT optimization of a real GeTe/TiTe 2 SL 23 resulting in a value of 3.4 Å for Te-Te interplanar distance and of 2.74 Å for Ge-Te interplanar distance.
As discussed in several previous works 21,24,25 , Te atoms are mostly three-fold coordinated in a pyramidal geometry, most of the Ge atoms are three-fold coordinated in a pyramidal geometry and four-or five-fold coordinated in a defective octahedral environment (octahedral bonding angles but coordination lower than six), while a minority fraction of Ge atoms are in tetrahedral geometries.The tetrahedral configuration is favored by homopolar Ge-Ge bonds 26 .These structural features are revealed by the bond angle distribution function (Fig. 2c): the peak at about 90 o in the angle distribution function of Ge atoms is due to pyramidal and defective octahedral configurations, the weak peak at about 170 o is due to axial bonds in a defective octahedral configuration and the shoulder at about 109 o is due to tetrahedra.A quantitative measure of the fraction of tetrahedral environments can be obtained from the local order parameter q introduced in Ref. 27 .It is defined as , where the sum runs over the pairs of atoms bonded to a central atom j and forming a bonding angle θ ijk .The order parameter evaluates to q = 1 for the ideal tetrahedral geometry, to q = 0 for the 6-fold coordinated octahedral site, to q = 5/8 for a 4-fold coordinated defective octahedral site, and q = 7/8 for a pyramidal geometry (three bonding angles of 90 • ).The distribution of the q order parameters for the slab and bulk models are compared in Fig. 2d.As discussed in previous works [28][29][30] , the fraction of Ge atoms in tetrahedral environments can be quantified by counting the 4-fold coordinated atoms with q in the range 0.8 < q < 1 .The resulting fraction of Ge atoms in tetrahedral geometry (with respect to the total number of Ge atoms) is 25 % in the slab and 19 % in the bulk.The coordination numbers of bulk a-GeTe at the experimental density of the crystalline α phase are larger than those of a-GeTe at the experimental density (0.03327 atoms/ Å 3 ) as expected due to the density increase (see for instance Ref. 21 ).
The amorphous model was then heated at six different target temperatures of 500, 550, 600, 650, 700, and 750 K to study the crystallization process, which was monitored in constant volume simulations, 1 ns long, at each temperature.The number of crystalline atoms was identified by the local order parameter Q dot 4 (i) (see Methods).Crystal nucleation starts at both surfaces of the amorphous slab, albeit the capping layers do not act as nucleation centers since there is no registry with the lattice of GeTe.
Snapshots of the crystallization process at 600 K are shown in Fig. 3; similar snapshots for the simulations at 650 and 700 K are given in Fig. S2 and Fig. S3 in the Supplementary Information.In all simulations we observed the formation of several crystalline nuclei.At the lower temperatures of 550-650 K, the final configuration is a polycrystalline model with several grains.At 700 and 750 K we observed the coarsening of different grains leading to a final single crystal model.The coarsening of the nuclei with different crystallographic orientation is a much slower process than the nucleation and growth of individual nuclei.Therefore, we observed coarsening only at the higher temperatures (700 and 750 K) because there the number of different nuclei is smaller due to the lower nucleation rate (the size of the critical nucleus is larger) and because the coarsening process  As already mentioned, the nucleation started at the surfaces, possibly because of a sort of layering due to the excess of Te on the outermost planes.Because of that, in the final configuration at 700 and 750 K where a single grain is formed, a Te layer is present at both surfaces which leads to the formation of 19 atomic layers instead of the 18 layers of the starting crystalline GeTe slab.As a consequence, several antisite defects (because of the Te plane is excess) and vacancies are present in the crystallized model.At all temperatures in the range 550-750 K, the crystalline phase is cubic-like as the β-phase of GeTe.By quenching to 300 K, the crystalline grains still keep a disordered β-phase configuration possibly because of a large number of defects.This occurs for all the models crystallized at different temperatures.
The number of crystalline atoms as a function of time for the simulations at different temperatures is shown in Fig. S4 in the Supplementary Information.The crystal growth velocity v g has been computed from the time derivative of the crystalline volume V c according to the scheme proposed in Ref. 32 as v g (t) = S −1 ac dV c /dt where S ac is the area of the crystal-amorphous interface (see Methods).The crystal growth velocity was calculated at different temperatures, after equilibration at the target temperature and before the different nuclei start to coalesce.The instantaneous v g , V c and S ac (see Ref. 31 ) as a function of time at the different temperatures is shown in Fig. 4. The average crystal growth velocities are obtained by averaging the instantaneous crystal growth velocity in a time interval of a few hundreds of ps as shown in Fig. S5 in the Supplementary Information.The resulting v g are compared in Table 2 with those of reference calculations for the bulk.We considered different models for the bulk amorphous phase at the density of crystalline α-GeTe.We first considered a slab of amorphous GeTe in contact with a slab of crystalline GeTe.The model was prepared starting from an orthorhombic cell of crystalline α-GeTe at the crystalline equilibrium density (0.0372 atoms/Å 3 ) consisting of 54 atomic layers with 180 atoms per layer.The system was then brought to 1500 K for 200 ps by letting 28 layers free to move and the other 26 frozen (model Bulk28 in Table 2).The molten region was then equilibrated at 1000 K for 100 ps and then quenched at 300 K in 100 ps, then all the atoms have set free to move and the system was equilibrated at the different target temperatures to compute the crystal growth velocities along the growth direction corresponding to the c axis of α-GeTe.
Moreover, we considered also a cubic 4096-atom bulk model of the amorphous phase still at the density of crystalline α-GeTe to study homogeneous crystal nucleation and growth.The number of crystalline atoms as www.nature.com/scientificreports/ a function of time at different temperatures for the two bulk models of a-GeTe are shown in Fig. S4 in in the Supplementary information.In the range 500-600 K, we see the formation of several overcritical nuclei, while at 650 K a single single nucleus forms and no nucleation occurs at and above 700 K up to 2 ns.The instantaneous v g , V c and S ac (see Ref. 31 ) as a function of time at the different temperatures is shown for the two bulk models in Fig. S6 and S7 in the Supplementary information.The average crystal growth velocities v g for the two bulk amorphous models are reported in Table 2.The results on v g for the bulk-like models are similar and about a factor two larger than those of the slab confined by the frozen TiTe 2 -like capping layers at the temperatures at maximal crystallization speed.The difference between the v g in the slab and in the bulk decreases by increasing temperature.The lower v g in the slab can be ascribed to the interaction among several nuclei which is enhanced with respect to the bulk because in the slab the nucleation centers appear just at the two surfaces.This means that the different nucleation centers are on average closer in the slab than in the bulk.There is a small reduction of the atomic mobility in the slab with respect to the bulk (see Table S1 and Fig. S8 in the Supplementary Information) which is, however, not sufficiently large to account for the difference in v g between the bulk and the slab.
Note that the crystal growth velocities of Table 2 for the bulk models are lower than those reported in previous simulations with the same NN potential because here the density is set to that of the α-GeTe which is higher than the experimental density of the amorphous phase considered in previous works 12 .A dependence of the crystal growth velocity on the density was also highlighted in previous works 33,34 (see for instance Fig. 2a in Ref. 33 ).
In the crystal growth from the amorphous/crystal interface in the bulk (Fig. S6 in the Supplementary Information), one observes an oscillation in time of the instantaneous crystal growth velocities.This can be ascribed to a different sticking coefficient of the Te and Ge atoms at the crystalline surface.In fact, the growth along the c axis of α-GeTe corresponds to the formation of the alternating layers made of only one atomic specie.The oscil- lation periodic is shorter than the time needed to complete a single crystalline layer because it results from the superposition of two growing surfaces.The same effect is also present at intermediate temperatures in the slab model (see Fig. S5 in the Supplementary Information).
Overall, we can conclude that the confinement has not a dramatic effect on the crystallization kinetic of our models of GeTe/TiTe 2 -like superlattice.Although a reduction in the crystal growth velocity with respect to the bulk is indeed observed at the lower temperatures, the maximum crystallization speed of the confined GeTe slab is comparable to that of the bulk at the same density which makes the GeTe/TiTe 2 superlattice a viable candidate for applications in neuromorphic computing with a foreseen improvement in data retention with respect to the Sb 2 Te 3 /TiTe 2 SL proposed in Ref. 7 .ac as described in Refs. 31.The v g reported in Table 2 are obtained by averaging the instantaneous v g over the time intervals highlighted in Fig. S5 in the Supplementary Information.Table 2. Crystal growth velocities v g (m/s) computed as described in Methods for the GeTe slab confined by the TiTe 2 -like capping layers (GeTe/TiTe 2 -like SL) for the bulk models (see text) with a thick (28 layers, Bulk28) molten region in contact with a crystalline slab and from homogeneous crystal nucleation and growth (Bulk homo ).is the reduced drift of the electrical resistance in the reset state due to reduced structural relaxations, i.e. aging, of the amorphous phase in confined ultrathin (about 5 nm thick) films.The same effect was also reported in Ref. 9 for ultrathin films of Sb (3-10 nm thick) confined by capping layers.It is therefore interesting to address whether a similar behavior could be observed in the GeTe/TiTe 2 superlattices.The drift in the electrical resistance is typically monitored on time scales (from seconds to hours) 7,35 not accessible by MD simulations.Nevertheless, we can attempt to estimate the effect of confinement on the structural relaxations in our models of the amorphous phase which are very far from the ideal glass because they are generated by fast quenching from the melt (100 ps).To this aim, starting from the models of bulk a-GeTe (4096-atom cubic box) and of the confined amorphous slab (about 3 nm thick) in the GeTe/TiTe 2 -like superlattice equilibrated at 300 K, we raised the temperature to 350 or 400 K in NVT simulations to accelerate the structural relaxations by still keeping the system below the glass transition to prevent crystal nucleation and growth.We considered models of bulk a-GeTe both at the experimental density of a-GeTe and the density of crystalline α-GeTe (see above).The energy gain due to the structural relaxations is then monitored as a function of time up to 4 ns as shown in Fig. 5. Actually, we do not observe sizeable differences in the behavior of a-GeTe in the bulk and in the confined slab which would suggest that the drift in a-GeTe is not mitigated by confinement in the geometry chosen here.The slightly faster structural relaxation in the sequence of bulk at the experimental density, GeTe/TiTe 2 SL and bulk at the crystal density is consistent with an acceleration of the drift by decreasing the density of the amorphous phase.We must, however, emphasizes once more that the time scale considered here for the structural relaxations (4 ns) is much shorter than the time scale over which the drift in the electrical resistance is measured experimentally for instance in the Sb 2 Te 3 /TiTe 2 superlattice of Ref. 7 (from seconds to hours), albeit the experimental data refer to 300 K while our simulations are performed at 350-400 K.

Conclusions
In summary, we performed MD simulations on the crystallization kinetics of ultrathin (3 nm) slabs of GeTe confined by capping layers aimed at mimicking the TiTe 2 spacers in so far hypothetical GeTe/TiTe 2 superlattices.This system is analog to the Sb 2 Te 3 / TiTe 2 heterostructure proposed in Ref. 7 for applications in neuromorphic devices.The replacement of Sb 2 Te 3 by GeTe would raise the crystallization temperature of the amorphous phase to improve the data retention.Since nanoconfinement is known to slow down the crystallization kinetics, for instance in elemental Sb 9 , we have investigated here nanoconfined GeTe to explore its potential applications in GeTe/TiTe 2 superlattices.The simulations show that nanoconfinement leads to a decrease of the crystal growth velocity v g with respect to the bulk amorphous phase at the same density, set here to that of the crystalline α phase.The reduction in v g , about a factor two at the temperature of maximal crystallization speed, is however rather minor in the perspective application in neuromorphic devices.As opposed to what found experimentally for Sb 2 Te 3 and elemental Sb 7,9 , we do not observe a reduction of structural relaxations, i.e. aging, in nanoconfined GeTe with respect to the bulk.However, we must remark that we have investigated structural relaxation in the amorphous phase on a time scale (4 ns at 400-350 K) much shorter than that considered experimentally to monitor the resistance drift (from seconds to hours at 300 K) in Sb 2 Te 3 /TiTe 2 superlattices 7 and in ultrathin GeTe layers (3 nm) in memory cells 35 .Therefore, we can not exclude that a mitigation of the drift could also be achieved in nanoconfined GeTe on the time scale of interest for the operation of the devices.In conclusions, MD simulation support the idea of investigating GeTe/TiTe 2 superlattices for applications in neuromorphic devices with improved data retention that we here put forward to future experimental work.q 4m,i q * 4m,j with q 4m (i) =

Figure 1 .
Figure 1.(a) Crystalline and (b) amorphous phase of the slab made of nine bilayers of α-GeTe encapsulated by capping layers.The capping layer is made by a frozen bilayer of crystalline GeTe at the lattice constant of TiTe 2 , aiming at mimicking the confining slabs of TiTe 2 in GeTe/TiTe 2 superlattices.Color code for atomic spheres: Ge (gray), Te (orange), Ge atoms in the capping layers (red).

Figure 2 .
Figure 2. Structural properties at 300 K of the slab of amorphous GeTe (black lines, see Fig.1b) confined by the TiTe 2 -like capping layers compared to those of a bulk model of amorphous GeTe at the density of α -GeTe generated by quenching from the melt (red lines).(a) Partial pair correlation functions.(b) Distribution of coordination numbers resolved per chemical species obtained by integrating the partial pair correlation functions up to the bonding cutoff of 3.00 Å for Ge-Ge and Te-Te bonds and 3.22 Å for Ge-Te bonds.(c) Bond angle distribution function resolved per central atomic species.The data are normalized to the number of triplets in each model.(d) Distribution of the q order parameter for tetrahedricity of the fourfold coordinated Ge atoms.

Figure 3 .
Figure 3. Simulation of the crystallization of a 3240-atom slab of amorphous GeTe capped by bilayers mimicking confinement by TiTe 2 in GeTe/TiTe 2 superlattices.Snapshots at different times at 600 K are shown for (a) 0.25 ns (b) 0.5 ns (c) 0.75 ns and (d) 1 ns.Crystallization starts at the surfaces of the amorphous slab, albeit the capping layers do not act as nucleation centers.Only crystalline atoms, identified by the Q dot 4 order parameter (see Methods), are shown.Different crystalline nuclei have different colors.

Figure 4 .
Figure 4. (a) Instantaneous crystal growth velocity v g , (b) volume occupied by the crystalline atoms V c and (c) area of the crystal-amorphous interface S ac as a function of time at the different temperatures in the crystallization of the superlattice configuration (GeTe/TiTe 2 -like SL).The crystal growth velocity is computed as v g = dV c /dtS −1 ac as described in Refs.31 .The v g reported in Table2are obtained by averaging the instantaneous v g over the time intervals highlighted in Fig.S5in the Supplementary Information.

Figure 5 .
Figure 5. Energy gain per atom as a function of time due to structural relaxations in NVT simulations at (a) 350 K and (b) 400 K for a-GeTe in a bulk model at the crystal density (red line), a-GeTe in a bulk model at the experimental density (orange line) and in confined slab geometry of the GeTe/TiTe 2 -like superlattice (green line).

Table 1 .
Average coordination number for different pairs of atoms computed from the partial pair correlation functions for the amorphous slab confined by the capping layers, compared with the data of a bulk amorphous model at the density of crystalline α-GeTe (see text).
Another advantage of the confinement in a SL geometry reported for Sb 2 Te 3 /TiTe 2 in Ref.